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An  improved  electrochemical  model  is  developed  to  study  the  ammonia  fed  solid  oxide  fuel  cell  based  on 
proton  conducting  electrolyte  (SOFC-H).  Including  the  chemical  reaction  kinetics  of  NH3  catalytic  thermal 
decomposition,  the  present  model  can  be  used  to  predict  the  performance  of  the  NH3  fed  SOFC-H  at 
an  intermediate  temperature  (i.e.  773  K).  Comparison  between  the  simulation  results  using  the  present 
model  and  experimental  data  from  literature  validates  the  accuracy  of  this  model.  Parametrical  analyses 
reveal  that  at  a  high  operating  temperature  (i.e.  1073  K),  the  NH3  fuel  is  completely  decomposed  to  H2  and 
N2  within  a  very  thin  layer  (30  (Jim)  near  the  anode  surface  of  an  SOFC-H.  It  is  also  found  that  operating 
the  NH3  fed  SOFC-H  at  an  intermediate  temperature  of  773  K  is  feasible  due  to  sufficiently  high  rate 
of  NH3  decomposition.  However,  further  decreasing  the  temperature  to  673  K  is  not  recommended  as 
less  than  10%  NH3  fuel  can  be  decomposed  to  H2  and  N2  in  the  SOFC-H.  The  effects  of  current  density 
and  electrode  microstructure  on  the  performance  of  the  NH3  fed  SOFC-H  are  also  studied.  It  is  found 
that  increasing  electrode  porosity  and  pore  size  is  beneficial  to  increase  the  partial  pressure  of  H2  at  the 
anode-electrolyte  interface.  The  model  developed  in  this  paper  can  be  extended  to  2D  or  3D  models  to 
study  practical  tubular  or  planar  SOFCs. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  have  been  identified  as  a  promis¬ 
ing  power  source  for  clean  and  efficient  electricity  generation  [1]. 
One  advantage  of  an  SOFC  is  its  fuel  flexibility.  As  SOFCs  work  at 
a  high  temperature  between  1073  K  and  1273  K,  various  alterna¬ 
tive  fuels,  such  as  hydrogen,  bio-methanol,  bio-ethanol  and  biogas 
can  be  directly  used  in  an  SOFC  [2-6].  In  recent  years,  ammonia 
(NH3)  emerges  as  a  promising  fuel  for  SOFCs  because  it  is  relatively 
cheap,  easy  to  store  and  transport,  and  relatively  safe  due  to  any 
leakage  being  easily  detectable  by  its  pungent  odor.  In  addition,  the 
infrastructure  of  ammonia  technology  has  been  well  established. 
Therefore,  there  is  increasing  interest  in  using  ammonia  in  fuel  cells, 
especially  SOFCs  [7-23]. 

Conventionally,  an  SOFC  employs  oxygen  ion  conducting  ceram¬ 
ics  as  its  electrolyte  (SOFC-O),  i.e.  yttria-stabilized  zirconia  (YSZ). 
In  order  to  achieve  high  performance,  SOFCs  are  usually  oper¬ 
ated  at  a  high  temperature  between  1073 1<  and  1273  K.  At  high 
temperatures,  both  electrochemical  reactions  and  ion  conduction 
are  fast,  which  in  turn  results  in  low  activation  and  ohmic  over¬ 


*  Corresponding  author.  Tel.:  +852  2859  2811;  fax:  +852  2858  5415. 
E-mail  addresses:  memni@graduate.hku.hk,  mengni@hku.hk  (M.  Ni). 

0378-7753 /$  -  see  front  matter  ©  2008  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2008.07.023 


potentials.  However,  the  high  operating  temperature  limits  the 
choice  of  materials  used  for  SOFCs  and  also  causes  catalyst  sin¬ 
tering  and  thermal  expansion  mismatch  of  SOFC  components.  In 
order  to  resolve  these  material  and  stability  problems,  a  lot  of 
research  efforts  have  been  made  to  operate  the  SOFC  at  intermedi¬ 
ate  temperatures.  However,  the  ionic  conductivity  of  YSZ  decreases 
considerably  with  decreasing  temperature,  which  can  cause  high 
ohmic  loss  at  the  electrolyte.  For  comparison,  some  proton  con¬ 
ducting  ceramics  show  good  ionic  conductivity  at  intermediate 
temperatures  (673-773  K)  and  are  thus  suitable  candidates  for  use 
as  SOFC  electrolyte.  Therefore,  there  is  increasing  interest  in  SOFC 
using  proton  conducting  electrolyte  (SOFC-H).  The  most  frequently 
used  proton  conducting  electrolytes  are  BaCe03 -based  ceramics, 
such  as  BaCe03  doped  with  Gd,  Nd,  Sm  or  Y  [14,15,21].  Compared 
with  conventional  SOFC-O,  the  SOFC-H  is  advantageous  as  complete 
fuel  utilization  can  be  easily  achieved  in  an  SOFC-H.  In  addition, 
when  ammonia  is  used  as  a  fuel,  an  SOFC-H  has  an  additional  advan¬ 
tage,  as  no  harmful  nitrogen  oxide  gas  will  be  generated  [7,8,20]. 
For  comparison,  nitrogen  oxide  may  be  formed  in  the  anode  cham¬ 
ber  of  an  SOFC-O  as  oxygen  ions  transporting  through  the  dense 
electrolyte  may  react  with  N2  at  the  surface  of  the  catalyst  particles. 
Previous  thermodynamic  analyses  have  demonstrated  that  fed  with 
hydrogen,  methane,  ethanol,  or  ammonia,  the  SOFC-H  has  a  higher 
maximum  energy  efficiency  than  the  SOFC-O  [24-28].  However, 
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the  above-mentioned  thermodynamic  studies  only  consider  the 
Nernst  potentials  of  SOFC-O  and  SOFC-H  while  the  overpotentials 
are  not  included.  In  a  recent  publication,  a  simplified  electrochemi¬ 
cal  model  has  been  developed  to  predict  the  actual  performance  of 
ammonia  fed  SOFC-H  [29].  The  model  considers  all  overpotentials 
involved  in  SOFC-H  operation  but  assumes  complete  decomposi¬ 
tion  of  NH3  at  the  anode  surface.  This  assumption  greatly  simplifies 
the  model  as  the  catalytic  thermal  decomposition  of  NH3  and  multi- 
component  mass  transfer  phenomena  in  the  porous  anode  can 
be  determined  separately.  However,  this  assumption  is  only  valid 
for  SOFC-H  working  at  a  high  temperature  (i.e.  1073 1<)  and  may 
become  invalid  for  SOFC-H  working  at  an  intermediate  tempera¬ 
ture  (i.e.  773 1<).  In  this  study,  an  improved  electrochemical  model 
is  developed  to  study  the  performance  of  the  NH3  fed  SOFC-H  at  an 
intermediate  temperature.  The  coupled  mass  transfer  and  chemical 
reaction  in  the  porous  anode  are  fully  considered  in  this  model. 


ature;  F  is  the  Faraday  constant;  PH2,a,  Ph2o,c  and  Po2,c  are  the 
partial  pressure  of  H2  (anode),  H20  (cathode),  and  02  (cathode), 
respectively.  The  superscript  ‘Int’  refers  to  the  electrode-electrolyte 
interface.  rjact<a  and  rjact<c  are  the  activation  overpotentials  at  the 
anode  and  cathode,  respectively,  which  can  be  calculated  with  the 
Butler-Volmer  equation  [31  ].  r]ohm  is  the  ohmic  overpotential  at  the 
electrolyte  and  can  be  analytically  calculated  with  Ohm’s  law  [31  ]. 
It  should  be  mentioned  that  the  concentration  overpotentials  are 
implicitly  included  in  the  reversible  potential  (F),  since  the  partial 
pressures  of  gaseous  species  used  in  Eq.  (3)  are  evaluated  at  the 
electrode-electrolyte  interface  [32].  In  order  to  obtain  these  inter¬ 
facial  partial  pressures,  the  characteristics  of  mass  transport  and 
chemical/electrochemical  reactions  in  the  porous  electrodes  must 
be  obtained. 

2.3.  Reaction  and  mass  transport  in  porous  electrode 


2.  Electrochemical  model 


2.  t.  Working  mechanisms  of  the  NH3  fed  SOFC-H 

The  working  mechanisms  of  the  NH3  fed  SOFC-H  are  shown  in 
Fig.  1.  At  the  anode,  NH3  undergoes  catalytic  thermal  decomposi¬ 
tion  over  a  Ni  catalyst  as 

2NH3  N2  +  3H2  (1) 

The  N2  is  removed  from  the  anode  chamber  by  the  gas  stream, 
while  the  H2  produced  undergoes  electrochemical  reactions  at 
the  anode-electrolyte  interface  to  produce  protons  and  electrons. 
The  electrons  traveling  from  the  anode  to  the  cathode  through  an 
external  circuit  yield  an  electrical  power  output.  The  protons  are 
transported  through  the  dense  electrolyte  to  the  cathode  chamber 
and  undergo  electrochemical  reactions  with  oxygen  molecules  and 
electrons  to  form  H20  at  the  cathode-electrolyte  interface. 


2.2.  Working  potential  of  the  NH3  fed  SOFC-H 


Considering  all  the  overpotentials  in  the  NH3  fed  SOFC-H,  the 
external  voltage  (V)  can  be  expressed  as 


V  —  E  —  ?7act,a  —  ?7act,c  —  Pohm 


F  —  Fq  + 
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where  Eq  is  the  ideal  standard  potential  of  the  SOFC-H,  which  can 
be  obtained  from  literature  [30,31  ];  E  is  the  reversible  (equilibrium) 
potential,  R  is  the  ideal  gas  constant;  T  is  the  operating  temper- 
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Fig.  1.  Working  mechanisms  of  the  NH3  fed  SOFC-H. 


2.3.1.  Anode 

In  literature,  there  are  a  lot  of  research  works  on  thermal  decom¬ 
position  of  ammonia  with  different  catalysts.  Many  metals  and 
compounds,  such  as  Fe,  Ni,  Pt,  Ru,  Ir,  Pd,  Ni-Pt,  Ni/Ru,  etc.  have 
been  examined  for  ammonia  thermal  decomposition  [33-42].  In 
most  of  these  studies,  dilute  ammonia  gas  (at  ppm  level)  was  used 
for  investigation.  Possible  reaction  steps  are:  (1 )  adsorption  of  NH3 
on  catalyst  surface;  (2)  cleavage  of  N-H  bond  of  adsorbed  NH3; 
and  (3)  recombinative  desorption  of  nitrogen  atoms  [37].  Although 
there  is  no  consensus  on  the  rate  limiting  steps  of  NH3  decomposi¬ 
tion  on  different  catalyst  and  under  different  operating  conditions, 
it  was  normally  observed  that  the  rate  of  ammonia  thermal  decom¬ 
position  is  inhibited  by  hydrogen  partial  pressure  at  a  relatively  low 
temperature  (below  673  K)  and  low  partial  pressure  of  NH3  (ppm 
level).  Based  on  experimental  observation,  Bradford  et  al.  [33  ]  found 
that  the  rate  of  NH3  decomposition  followed  the  equation: 

r  =  kpm3K2  (4) 

where  a  varied  from  0.69  to  0.75;  varied  from  -1.5  to  -2.0;  and  k 
is  a  constant  for  Eq.  (4).  The  reported  activation  energy  was  about 
23kcalmol_1.  This  above  equation  reveals  the  strong  inhibition 
effect  of  hydrogen  partial  pressure  on  the  thermal  decomposition  of 
dilute  NH3  (ppm  level)  at  a  temperature  of  below  673  K.  In  a  recent 
study  of  NH3  thermal  decomposition  at  a  temperature  of  between 
623 1<  and  923 1<  and  a  NH3  partial  pressure  between  0.5  kPa  and 
2.0  kPa,  the  reaction  rate  of  NH3  decomposition  over  Ru  based  cat¬ 
alyst  was  found  to  satisfy  the  following  equation  [34]: 

r  =  ^nh3  (5) 

where  a  was  about  2.0,  indicating  that  the  reaction  rate  was  sec¬ 
ond  order  with  respect  to  NH3  partial  pressure.  The  second  order 
dependence  on  NH3  partial  pressure  was  explained  by  the  assump¬ 
tion  that  the  recombinative  desorption  of  nitrogen  atoms  was  the 
rate-determining  step  [34,35]. 

The  above  results  are  applicable  to  dilute  NH3  at  a  ppm  level. 
Chellappa  et  al.  [36]  investigated  the  decomposition  kinetics  of  pure 
ammonia  over  Ni-Pt/Al203.  They  found  that  the  decomposition 
kinetics  of  pure  ammonia  was  quite  different  from  those  reported 
previously  with  diluted  ammonia  at  a  ppm  level.  At  a  temperature 
between  793  K  and  963  K,  the  reaction  rate  was  found  first  order 
with  respect  to  NH3  partial  pressure  (Eq.  (6))  and  no  H2  inhibi¬ 
tion  effect  was  observed.  In  addition,  the  activation  energy  (Fa)  was 
found  to  be  about  46.9  kcal  mol-1  (equivalent  to  196.2  kj  mol-1 ). 

r  =  k0exp(-||)pNH3  (6) 

Obviously,  the  concentration  of  NH3  in  the  NH3  fed  fuel  cells 
should  be  high  (much  higher  than  the  ppm  level),  thus  the  first 
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order  dependence  of  reaction  rate  on  partial  pressure  of  NH3  was 
adopted  in  this  research.  The  inhibition  effect  of  H2  partial  pressure 
is  assumed  to  be  negligible  as  it  is  significant  only  at  low  tempera¬ 
ture  and  very  low  NH3  concentration  [37].  Recently,  Liu  et  al.  [38] 
investigated  the  decomposition  kinetics  of  NH3  over  nano-sized 
Ni/SBA-15  (SBA:  Santa  Barbara  Amorphous)  catalysts.  The  rate  of 
H2  formation  (mmol  min-1  g“a\)  increased  from  about  8.0  at  723  K 
to  about  33  at  923  K.  At  a  typical  temperature  of 773  K,  the  H2  forma¬ 
tion  rate  was  about  17.5  mmol  min-1  g-^.  The  pore  volume  of  the 
Ni/SBA  based  catalysts  ranged  from  0.36  cm3  g-1  to  0.82  cm3  g-1. 
In  another  recent  study  conducted  by  Zheng  et  al.  [39],  the  kinet¬ 
ics  of  NH3  decomposition  to  H2  over  Ni/Al203  catalysts  have  been 
investigated.  It  was  reported  that  at  a  temperature  of  773  K,  the  rate 
of  H2  formation  was  about  24.1  mmol  min-1  g-^.  The  pore  volume 
of  the  Ni/Al203  based  catalysts  was  found  to  be  about  0.3  cm3  g-1. 
Assuming  a  porosity  of  30%,  the  volume  of  the  porous  catalyst  was 
about  1.0  cm3  g-1 .  In  addition,  other  investigators  also  reported  that 
the  rate  of  hydrogen  formation  was  about  20  mmol  min-1  g-^  at  a 
temperature  of  773  K  [40-42].  Adopting  this  hydrogen  formation 
rate  at  773  K,  the  coefficient  k0  can  thus  be  determined  (4.0  x  1015) 
from  Eq.  (6).  The  rate  of  NH3  decomposition  can  thus  be  expressed 
as 


Df 


0.00133, 


(±  JL 

+  Mj 


0.5 


jl.5 


PafjG  d 


(14) 


where  £,  £,  and  rp  are  the  porosity,  tortuosity  and  mean  pore  radius 
of  the  electrodes;  Mj  is  the  molecular  weight  of  species  i ;  C2D  is 
a  dimensionless  diffusion  collision  integral;  and  o jj  is  the  mean 
characteristic  length  of  species  i  and  j.  The  values  of  and  o jj  can 
be  obtained  from  the  literature  [48]. 

The  pressure  gradient  (dP/dx)  can  be  evaluated  with  a  method 
developed  by  Zhu  and  Kee  [49]; 


dp  ELiW/pff) 

dx  ( 1  /RT)  +  (BoP/RTm)]^  (yt/DJ$) 


At  the  anode  surface,  the  molar  fractions  of  each  species 
are  known.  At  the  electrode-electrolyte  interface,  electrochemical 
reactions  take  place  and  thus  the  flux  of  H2  can  be  related  with  the 
current  density  as 


Nh2 
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r  =  4.0x10- exp  (-1^2)  PNH3  (7) 

In  a  steady  state,  the  transport  of  each  participating  component 
is  determined  by  the  local  conservation  of  mass: 


where  J  is  the  current  density  (Am-2)  and  da  is  the  thickness  of 
the  anode.  As  the  remaining  species  are  not  involved  in  the  electro¬ 
chemical  reaction,  their  fluxes  at  the  electrode-electrolyte  interface 
are  equal  to  zero. 


dNNH3 
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where  Nj  is  the  flux  of  species  i  (mol  m-2  s-1 );  x  is  the  depth  mea¬ 
sured  from  the  electrode  surface.  The  subscript  ‘a’  refers  to  the 
anode. 

The  dusty  gas  model  (DGM)  is  used  to  model  the  multi- 
component  mass  transfer  within  the  anode  and  cathode  as  it  can 
predict  the  mass  transfer  in  porous  media  more  accurately  than 
the  Fields  model  and  Stefan-Maxwell  model  [43,44].  According  to 
DGM,  the  transport  of  gas  species  in  both  anode  and  cathode  is 
governed  by 
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where  Nj  is  the  flux  of  gas  species  i ;  y2-  is  the  molar  fraction  of 
gas  species  i;  Df  £  is  the  effective  Knudsen  diffusion  coefficient  of 
species  i;  D?.ff  is  the  effective  binary  molecular  diffusion  coefficient 
of  species  i  and  j;  fi  is  the  viscosity  of  the  gas  mixture;  and  B0  is 
the  permeability  of  the  porous  electrode,  which  can  be  calculated 
by  the  Kozeny-Carman  relationship  [44]: 


(NNH3Lda=°  3nd  NN2,a|x=da  =0) 

The  above  governing  equations  (Eqs.  (7)-(15))  are  inter-related 
differential  equations  and  can  be  solved  by  numerical  method. 
Finite  difference  method  was  used  to  discretize  the  governing  equa¬ 
tions.  An  iterative  scheme  was  developed  to  obtain  the  solution  at 
the  discrete  grids.  Calculations  were  repeated  till  convergence  was 
achieved.  Solving  the  above-mentioned  equations,  the  partial  pres¬ 
sure  of  H2  at  the  anode-electrolyte  interface  can  thus  be  obtained. 


2.3.2.  Cathode 

In  the  cathode,  02  is  electrochemically  consumed  while  H20  is 
electrochemically  produced  and  no  chemical  reaction  is  involved. 
Therefore,  the  transport  of  gas  species  (02,  N2  and  H20)  can  be 
determined  by 

hm. 

-j-I=0  (i  =  02,  N2  and  H20)  (17) 

The  fluxes  of  the  transporting  species  can  be  determined  by  the 
DGM  (Eq.  (11)).  Similar  to  the  anode,  the  molar  fractions  of  gas 
species  at  the  cathode  surface  are  given.  At  the  cathode-electrolyte 
interface,  the  fluxes  of  02  and  H20  can  be  related  with  the  current 
density  as 
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The  effective  diffusion  coefficients  (Df£,  D?.ff)  can  be  determined  by 
the  following  equations  [45-47]: 
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After  solving  the  governing  equation  (Eqs.  ( 11 )  and  ( 15))  with  the 
above  boundary  conditions,  the  partial  pressures  of  02  and  H20  at 
the  cathode-electrolyte  can  thus  be  obtained.  After  obtaining  the 
partial  pressures  of  electrochemical  reaction  species  (H2,  H20  and 
02 )  at  the  electrode-electrolyte  interface,  the  equilibrium  potential 
of  the  NH3  fed  SOFC-H  can  thus  be  calculated  with  Eq.  (3). 
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Table  1 

Input  parameters  for  the  present  electrochemical  model  of  the  NH3  fed  SOFC-H 


Parameter 

Value 

Operating  temperature,  T  (K) 

1073,  873,  773,  673 

Operating  pressure,  P  (atm) 

1.0 

Gas  composition  at  the  electrode  surface  (%) 

Anode:  NH3 

100 

Cathode:  H20/02/N2 

3/18/79 

Exchange  current  density  of  anode 

5300.0  [30,31] 

(Am-2) 

Exchange  current  density  of 

2000.0  [30,31] 

cathode  (Am-2) 

Electrode  porosity 

0.4  [45] 

Electrode  tortuosity 

5.0  [45] 

Electrode  pore  radius  (|Jim) 

0.5  [45] 

Conductivity  of  electrolyte  (i.e. 

0.94  Sm-1  (estimated  from  [50]) 

Sm-doped  BaCe03)  at  1073  K 

Component  thickness  (p,m) 

Electrolyte 

50.0 

Anode 

500.0 

Cathode 

50.0 

3.  Results  and  discussion 

3  A.  Model  evaluation 

In  the  literature,  there  are  very  few  experimental  data  avail¬ 
able  on  the  NH3  fed  SOFC-H.  In  this  section,  an  evaluation  of 
the  model  was  conducted  by  comparing  the  present  modeling 
results  with  the  experimental  data  from  literature  [8].  In  the  exper¬ 
iment,  the  Gd-doped  BaCe03  (BCGO)  electrolyte  was  prepared  by 
a  modified  glycine-nitrate  process  with  stoichiometric  precursors 
of  Ba(CH3COO)2,  Ce(N03)2  and  Gd(N03)2.  The  composite  anode 
substrate  and  the  BCGO  electrolyte  film  were  fabricated  by  the 
dry-pressing  method.  The  SOFC-H  was  based  on  anode-supported 
configuration,  in  which  the  thickness  of  anode  and  electrolyte 
were  650  p,m  and  50  pan,  respectively.  Commercial  industrial  grade 
ammonia  and  humid  oxygen  (3%  H20)  was  used  as  the  fuel  and 
oxidant,  respectively.  The  tests  were  conducted  at  a  typical  tem¬ 
perature  of  973  K  and  a  pressure  of  1.0  atm.  The  values  of  input 
parameters,  such  as  porosity,  pore  radius,  and  tortuosity,  are  sum¬ 
marized  in  Table  1 .  As  can  be  seen  from  Fig.  2,  the  present  simulation 
results  agreed  reasonably  well  with  the  experimental  data.  There¬ 
fore,  the  model  developed  in  this  paper  is  reliable  for  the  following 
parametrical  analyses. 


Fig.  2.  Comparison  between  simulation  results  with  experimental  data  from  liter¬ 
ature  for  model  validation. 


3.2.  Temperature  effect 

Temperature  is  an  important  parameter  governing  the  per¬ 
formance  of  the  NH3  fed  SOFC-H.  In  this  section,  the  effects  of 
temperature  on  the  NH3  catalytic  thermal  decomposition  kinetics 
and  the  SOFC-H  performance  are  investigated. 

Fig.  3  shows  the  temperature  dependence  of  NH3  catalytic  ther¬ 
mal  decomposition  at  a  NH3  pressure  of  1.0  atm.  As  can  be  seen, 
the  rate  of  NH3  thermal  decomposition  strongly  depends  on  tem¬ 
perature,  i.e.  the  reaction  rate  increases  significantly  from  about 
1.0  mol  m-3  s-1  at  a  temperature  of  673  K  to  3.6  x  107  mol  m-3  s-1 
at  a  temperature  of  1273  K.  This  is  consistent  with  the  experimental 
observations  in  literature  [36,37]. 

The  distributions  of  NH3  catalytic  thermal  decomposition  rate  in 
the  porous  anode  of  an  SOFC-H  at  various  temperatures  are  shown 
in  Fig.  4(a)-(d).  For  SOFC-H  operated  at  a  high  temperature  (1073  K 
or  873  K),  the  rate  of  NH3  catalytic  thermal  decomposition  is  high 
and  the  NH3  fuel  can  be  completely  decomposed  before  reaching 
the  anode-electrolyte  interface  (Fig.  4(a),  (b),  and  (e)).  Especially  at 
a  temperature  of  1073  K,  the  rate  of  NH3  decomposition  decreases 
very  quickly  from  about  4.0  x  105  molm-3  s-1  to  0.0  mol  m-3  s-1 
within  a  very  thin  layer  (about  30  |jim)  near  the  anode  surface 
(Fig.  4(a)).  This  is  because  the  rate  of  NH3  decomposition  is  so  high 
that  the  NH3  fuel  is  completely  decomposed  to  H2  and  N2  in  a  thin 
layer  of  the  anode  surface  (Fig.  4(e)).  In  an  experimental  investiga¬ 
tion  of  the  NH3  fed  SOFC-H  at  a  temperature  of  1073  K,  Zhang  et  al. 
observed  that  NH3  was  decomposed  to  H2  and  N2  before  reaching 
the  anode-electrolyte  interface  [18].  Therefore,  the  present  finding 
agrees  with  the  experimental  observations  in  literature.  In  our  pre¬ 
vious  electrochemical  model  of  the  NH3  fed  SOFC-H  working  at  a 
temperature  of  1073  K,  we  assume  that  NH3  is  completely  decom¬ 
posed  near  the  anode  surface  [29].  The  finding  of  this  study  justifies 
the  assumption  made  in  our  previous  study.  When  the  SOFC-H  is 
operated  at  a  lower  temperature  (i.e.  773  K  or  673 1<),  the  rate  of  NH3 
decomposition  is  lower  and  only  part  of  the  NH3  fuel  can  be  decom¬ 
posed  to  H2  and  N2  in  the  SOFC-H  (Fig.  4(c)-(e)).  At  a  temperature 
of  673  K,  less  than  10%  NH3  fuel  can  be  decomposed  throughout  the 
anode  of  SOFC-H  because  of  too  low  NH3  decomposition  kinetics 
(Fig.  4(d)  and  (e)).  From  experiments,  it  is  known  that  the  perfor¬ 
mance  of  the  NH3  fed  SOFC-H  is  highly  dependent  on  the  extent  of 
NH3  decomposition  in  the  SOFC-H  [8,9,11  ].  Therefore,  it  can  be  seen 
from  the  present  study  that  a  too  low  temperature  (i.e.  673 1<)  is  not 
feasible  for  the  NH3  fed  SOFC-H,  unless  very  active  catalyst  for  NH3 
decomposition  can  be  developed.  On  the  other  hand,  a  considerable 
amount  of  NH3  can  be  decomposed  to  H2  and  N2  at  a  temperature  of 
773  K  in  the  SOFC-H,  although  the  conversion  is  not  100%  (Fig.  4(c) 
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Fig.  4.  Effect  of  operating  temperature  on  the  performance  of  the  NH3  fed  SOFC-H:  (a)  NH3  decomposition  rate  at  1073  K;  (b)  NH3  decomposition  rate  at  873  K;  (c)  NH3 
decomposition  rate  at  773  K;  (d)  NH3  decomposition  rate  at  673  K;  and  (e)  molar  fraction  of  NH3. 


and  (e)).  Therefore,  it  is  possible  to  operate  the  NH3  fed  SOFC-H  at 
a  temperature  of  773  K  with  acceptable  cell  performance.  In  Ma’s 
study,  it  is  reported  that  100%  NH3  decomposition  is  achievable  for 
SOFC-H  working  at  a  temperature  of  873  K  while  the  decomposition 
of  NH3  is  incomplete  when  the  temperature  is  reduced  to  773  K  [9]. 
In  another  study  conducted  by  Maffei  et  al.  [17],  it  is  found  that  the 
NH3  conversion  is  close  to  100%  at  a  temperature  of  873  K,  while  it 
is  about  30-55%  at  a  temperature  of  773  K.  The  present  modeling 
results  are  consistent  with  experimental  data  in  literature. 

In  the  following  analyses,  the  temperature  of  773 1<  is  selected 
as  a  typical  intermediate-temperature  for  the  NH3  fed  SOFC-H. 
For  any  fuel  fed  SOFC-H,  the  porous  cathode  involves  the  same 
electrochemical  reaction  and  mass  transport  phenomena,  i.e.  02  is 
consumed  while  H20  is  electrochemically  produced  at  the  porous 


cathode.  As  the  SOFC-H  fed  with  H2  and  CH4  has  been  investi¬ 
gated  in  our  previous  studies  [45,51  ],  the  following  sections  will  be 
focused  on  the  chemical  reaction  and  mass  transport  in  the  porous 
anode  of  the  NH3  fed  SOFC-H. 

3.3.  Current  density  effect 

Fig.  5  shows  the  effect  of  current  density  on  the  NH3  decom¬ 
position  kinetics  and  the  SOFC-H  performance  at  the  typical 
temperature  of  773  K.  The  pressure  distribution  in  the  porous  anode 
of  the  NH3  fed  SOFC-H  is  shown  in  Fig.  5(a).  At  a  low  current 
density  (i.e.  2000 Am-2),  the  anode  pressure  initially  increases 
with  increasing  anode  depth,  and,  after  reaching  a  maximum, 
decreases  slightly  with  further  increasing  anode  depth  (Fig.  5(a)). 
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Fig.  5.  Current  density  effect:  (a)  pressure  in  the  anode;  (b)  molar  fraction  of  H2  in  the  anode;  (c)  rate  of  NH3  decomposition;  and  (d)  partial  pressure  of  H2  at  the 
anode-electrolyte  interface. 


The  molar  fraction  of  H2  shows  a  similar  dependence  on  the  anode 
depth  (Fig.  5(b)).  There  are  two  mechanisms  governing  the  pres¬ 
sure  in  the  anode:  (1)  decomposition  of  NH3  to  H2  and  N2,  which 
tends  to  increase  the  pressure  and  (2)  consumption  of  H2  at  the 
anode-electrolyte  interface  due  to  electrochemical  reactions  (H20 
is  electrochemically  produced  at  the  cathode),  which  tends  to 
reduce  the  pressure  in  the  anode.  At  the  surface  of  the  anode,  the 
decomposition  of  NH3  is  important  due  to  a  high  partial  pressure 
of  NH3.  Therefore,  the  pressure  increases  with  increasing  anode 
depth  near  the  anode  surface.  With  an  increase  in  anode  depth, 
the  rate  of  NH3  decreases  due  to  a  lower  partial  pressure  of  NH3 
(multiplication  of  the  anode  pressure  and  the  NH3  molar  frac¬ 
tion)  (Fig.  5(c)).  Thus,  the  pressure  tends  to  decrease  in  a  deeper 
layer  because  the  resistance  of  the  porous  structure  to  the  trans¬ 
port  of  gaseous  species  is  relatively  more  important  (Fig.  5(a)).  At 
higher  current  densities  (i.e.  5000  A  m-2  or  8000  A  m-2 ),  the  anode 
pressure  decreases  monotonically  with  increasing  anode  depth, 
indicating  that  the  transporting  resistance  is  the  dominating  mech¬ 
anisms  governing  the  anode  pressure.  Obviously,  the  molar  fraction 
of  H2  decreases  with  increasing  current  density  (Fig.  5(b)),  which 
in  turn  results  in  a  higher  molar  fraction  of  NH3.  The  combined 
effect  of  higher  NH3  molar  fraction  and  lower  anode  pressure  at 
higher  current  density  leads  to  a  slightly  higher  NH3  decomposi¬ 
tion  rate  at  a  higher  current  density  (Fig.  5(c)).  In  quantifying  the 
electrochemical  performance  of  an  SOFC,  the  partial  pressure  of  H2 
at  the  anode-electrolyte  interface  is  often  an  important  parame¬ 
ter,  as  it  is  directly  related  to  the  working  potential  of  an  SOFC. 
As  both  the  anode  pressure  and  the  H2  molar  fraction  decrease 
with  increasing  current  density,  the  partial  pressure  of  H2  at  the 
anode-electrolyte  interface  is  found  to  decrease  with  increasing 
current  density  (Fig.  5(d)),  which  is  consistent  with  previous  anal¬ 
ysis  on  SOFC-O  and  SOFC-H  [29,31,45]. 


3.4.  Micro- structural  effect 

In  addition  to  the  above-mentioned  operating  parameters  (tem¬ 
perature  and  current  density),  the  electrode  microstructures  also 
significantly  influence  the  performance  of  SOFC-H.  In  this  section, 
the  effect  of  electrode  microstructures  on  the  NH3  fed  SOFC-H  is 
investigated  at  a  typical  current  density  of  5000  Am-2  and  a  typ¬ 
ical  temperature  of  773  K.  Electrode  porosity  and  pore  radius  are 
studied  due  to  their  influences  in  SOFC  performance. 

Fig.  6(a)  shows  the  distribution  of  pressure  in  the  porous 
anode  with  varying  electrode  porosities.  The  anode  pressure  is 
found  to  increase  with  an  increase  in  electrode  porosity.  This  is 
because  the  diffusion  coefficients  (Eqs.  (13)  and  (14))  are  higher 
at  a  higher  electrode  porosity.  In  addition,  the  H2  molar  fraction 
is  slightly  higher  at  a  higher  electrode  porosity,  which  in  turn 
causes  the  NH3  molar  fraction  to  decrease  slightly  with  increasing 
electrode  porosity.  As  the  anode  pressure  increases  considerably 
with  increasing  electrode  porosity  and  the  reduction  in  NH3  molar 
fraction  is  insignificant,  the  rate  of  NH3  decomposition  is  found 
higher  at  a  higher  electrode  porosity  (Fig.  6(b)).  Due  to  the  slight 
increase  in  H2  molar  fraction  and  a  considerable  increase  in  the 
anode  pressure,  the  partial  pressure  of  H2  at  the  anode-electrolyte 
interface  increases  with  increasing  electrode  porosity  (Fig.  6(c)). 
The  effect  of  electrode  pore  radius  on  the  anode  pressure  of  the 
NH3  fed  SOFC-H  is  shown  in  Fig.  6(d).  Similar  to  the  electrode 
porosity,  the  anode  pressure  increases  with  increasing  electrode 
pore  radius  due  to  higher  diffusion  coefficients  (Eqs.  (13)  and 
(14))  and  a  lower  transport  resistance.  The  rate  of  NH3  decom¬ 
position  also  increases  with  increasing  electrode  pore  radius  due 
to  enhanced  partial  pressure  of  NH3  (Fig.  6(e)).  In  addition,  the 
partial  pressure  of  H2  at  the  anode-electrolyte  interface  initially 
increases  considerably  with  increasing  electrode  pore  radius  and 
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Depth  in  anode  /  p/m  Electrode  pore  radius  /  pm 


Fig.  6.  Micro-structural  effect:  (a)  effect  of  electrode  porosity  on  anode  pressure;  (b)  effect  of  electrode  porosity  on  NH3  decomposition  rate;  (c)  effect  of  electrode  porosity 
on  partial  pressure  of  H2  at  the  anode-electrolyte  interface;  (d)  effect  of  electrode  pore  radius  on  anode  pressure;  (e)  effect  of  electrode  pore  radius  on  NH3  decomposition 
rate;  and  (f)  effect  of  electrode  pore  radius  on  partial  pressure  of  H2  at  the  anode-electrolyte  interface. 


then  tends  to  level  off  when  the  electrode  pore  radius  is  higher 
than  0.6  \xm  (Fig.  6(f)).  This  is  because  at  an  electrode  pore  radius 
higher  than  0.6  pan,  both  H2  molar  fraction  and  the  anode  pres¬ 
sure  only  slightly  increase  with  a  further  increase  in  electrode  pore 
radius. 

4.  Conclusion 

An  improved  electrochemical  model  was  developed  to  study  the 
performance  of  the  NH3  fed  SOFC-H.  This  model  is  an  extension 
of  our  previous  model  assuming  complete  decomposition  of  NH3 
at  the  anode  surface.  Including  the  kinetics  of  NH3  catalytic  ther¬ 
mal  decomposition,  the  present  model  is  capable  of  predicting  the 
performance  of  the  NH3  fed  SOFC-H  at  intermediate  temperatures 


(i.e.  773  K).  The  simulation  results  are  compared  with  experimental 
data  from  literature  and  good  agreement  is  found  between  them, 
thus  the  model  presented  in  this  paper  is  validated. 

Parametrical  simulations  are  performed  to  study  in  detail  the 
performance  of  the  NH3  fed  SOFC-H.  It  is  found  that  the  NH3 
catalytic  decomposition  is  highly  temperature-dependent.  At  a 
temperature  of  1073  K,  the  NH3  fuel  can  be  completely  decom¬ 
posed  within  a  very  thin  layer  (30  pan)  near  the  anode  surface. 
With  a  decrease  in  operating  temperature,  the  rate  of  NH3  cat¬ 
alytic  thermal  decomposition  decreases  considerably.  It  is  found 
that  it  is  feasible  to  operate  the  NH3  fed  SOFC-H  at  an  inter¬ 
mediate  temperature  of  773  K,  however,  further  reduction  of  the 
temperature  to  673  K  will  dramatically  decrease  the  SOFC-H  per¬ 
formance  as  less  than  10%  NH3  can  be  decomposed  to  H2  and  N2. 
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Increasing  current  density  slightly  increases  the  rate  of  NH3  decom¬ 
position  but  considerably  decreases  the  partial  pressure  of  H2  at  the 
anode-electrolyte  interface.  In  addition,  increased  electrode  poros¬ 
ity  and  pore  size  are  both  found  beneficial  to  increase  the  partial 
pressure  of  H2  at  the  anode-electrolyte  interface. 

The  result  of  this  study  justifies  the  assumption  made  in  our 
previous  publication  assuming  complete  NH3  decomposition  at  the 
surface  of  SOFC-H  anode.  The  model  developed  in  this  paper  can  be 
used  to  predict  the  performance  of  the  NH3  fed  SOFC-H  for  design 
optimization. 
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